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Abstract Motivated by the wide applications of BTZ black holes and interesting results of gravity’s 
rainbow, we consider three dimensional rainbow solutions and investigate their thermodynamic prop¬ 
erties. In addition to investigate black holes thermodynamics related to AdS/CFT correspondence, one 
may regard gravity’s rainbow to to encode quantum gravity effects into the black hole solutions. We 
take into account the various models of linear and nonlinear electrodynamics and study their effects 
on the gravity’s rainbow spacetime. We also examine thermal stability and find that obtained three 
dimensional rainbow black holes are thermally stable. 


1 Introduction 

One of the weighty dream of physicists is creating a consistent quantum theory of gravity. Although 
various aspects of the quantum gravity is still very much alive, till now, all the attempts of obtaining 
a complete description of the quantum gravity have been unfeasible. Gravity’s rainbow arises from 
deep attempts to fix the gaps between the theory of general relativity and quantum mechanics. One 
of the possible interesting consequence of gravity’s rainbow is that the universe had no beginning, and 
therefore time is stretched back infinitely without encountering a big bang [I]. Also, thermodynamical 
properties of black p-branes was reported within the framework of gravity’s rainbow [2J. Furthermore, 
as more recent application of gravity’s rainbow, a recent investigation on the possibility of resolving 
the black hole information paradox was carried out in Ref. 0]. In addition, the effects of the gravity’s 
rainbow on the deflection of light, photon time delay, gravitational red-shift, and the weak equivalence 
principle have been studied in Ref. [1]. Other consequences of the gravity’s rainbow in the context of 
the gravitational collapse and black hole physics have been investigated in literature (see IMEIGIIB] and 
references therein). 

On the other hand, some of theoreticians believe that the violation of Lorentz invariancy is an es¬ 
sential rule to construct quantum theory of gravity. The Lorentz invariance violation may be expressed 
in form of modified dispersion relations imnnumiimiis]- In other words, according to the loop quantum 
gravity (LQG) results JTO] and spacetime discreteness assumption [13] at the Planck scale, one may 
regard Lorentz invariance violation or modified dispersion relation by redefining the physical momen¬ 
tum and energy. In order to have a deep insight into the Lorentz invariancy, one may regard the close 
relation between gravity’s rainbow and Horava-Lifshitz gravity M- This relation comes from the fact 
that both these theories are based on the breaking of the usual energy-momentum dispersion relation 
in the UV limit in such a way that it reduces to the usual energy-momentum dispersion relation in 
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the IR limit. The UV modification of the usual energy-momentum relation implies the breaking of the 
Lorentz symmetry. Due to the existence of an unstable perturbative string vacuum, the spontaneous 
breaking of the Lorentz symmetry can be occurred in string theory [51. In addition, the Lorentz sym¬ 
metry breaking can be studied using black brane in the context of type IIB string theory m- Taking 
into account the mentioned motivation of UV deformation of geometries that occur in string theory 
and the close relation between Horava-Lifshitz gravity and gravity’s rainbow m > recently rainbow de¬ 
formation of geometries has been performed, in which leads to modification of the dispersion relation. 
One of the modified dispersion relation, which is predicted by new approach of quantum gravity mi 
ezi, takes the following general form 


f 2 (s)E 2 -g 2 (e)p 2 =m 2 , (1) 

where dimensionless relative energy e = E/E p , in which E is the energy of the probing particle and 
E p is the Planck energy. In addition, temporal and spatial rainbow functions (/(e) and g(e)) satisfy 
the following conditions 

hm/(e) = l, limo(e) = 1. (2) 

e-s-0 e->0 

On the other hand, such modified dispersion relation can be originated from the so-called doubly 
special relativity 1181 . Doubly special relativity is a natural extension of special relativity, which enjoys 
the invariancy of the speed of light, to the case of assuming the invariancy of the maximum energy scale 
(the Planck energy) or equivalently minimum length scale [151 . Besides, the generalized uncertainty 
principle (GUP) which comes from various approaches to quantum gravity, is based on the existence of 
minimum measurable length scale pnil2 l 22 23.24. 2 51I3H] . Such minimum length scale has been applied 
to black hole thermodynamics [27ll281f29ll301131] . From string theory point of view, it is not possible to 
probe spacetime below the string length scale and one can translate the minimum length scale into a 
maximum energy scale Including the curvature into the doubly special relativity, 

Magueijo and Smolin proposed the so-called doubly general relativity [31]. In other words, doubly 
general relativity is a natural extension of doubly special relativity to the case of assuming curved 
spacetime. Regarding the doubly general relativity, one finds the geometry of spacetime may depend 
on the energy of the particle moving in it. In other words, spacetime is parameterized by the ratio 
e =E/Ep to obtain a parametric family of metrics, the so-called rainbow of metrics [38]. Hence, the 
modified metric in gravity’s rainbow can be written as l38| 

g(e) = g ab e a {e) ® e b {e), (3) 


in which the energy dependence of the frame fields is 


eo(e) = 




(4) 


where the hatted quantities refer to the energy independent frame fields. In Ref. [8] it was shown 
that one can explain the absence of black holes at LHC by using the gravity’s rainbow. In addition, 
regarding gravity’s rainbow, one may resolve the black hole information paradox without any need for 
black hole complementarity or Firewall [3]. The functional forms of f(e) and g{e) are based on various 
phenomenological motivations mmmmmmmm ■ Among different proposals in the literature, 
we point out three interesting models of modified dispersion relation. 

One of the most studied models was proposed in Refs. mm- This model, which we identify it as 
first model, is compatible with the results of LQG and non-commutative spacetime m ■ In this model 
one should set the temporal rainbow’s function to unity and spatial one is a square root function with 
the following explicit forms 

/ (e) = 1, 9 (e) = V 1 - V £n ) (5) 

Amclino-Camelia, et al. proposed another model in Ref. [T2], which we will call second model. The 
hard spectra of gamma-ray bursters may be explained by the second model ]12] ■ In this model one sets 
spatial rainbow function to unity and temporal one has exponential form 


/(e) = 


e' 76 - 1 


( 7 £ 


9 (e) = 1, 


( 6 ) 
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Another model, which we will recognize it as third model, is related to constant speed of light and 
one may use it to solve the horizon problem |3S.. In this model both rainbow’s functions are the same 
with the following fractional form 

/(e)=s,£) = rb' (7) 

Regarding the three mentioned models, it is notable that the dimensionless parameters rj , n, a and 
A are of the order of unity. 

Among all black hole solutions one can think of, BTZ (Banados- Teitelboim-Zanelli) H5H9] black 
hole can be considered as the unique one. It is believed that (2 + 1)—dimensional BTZ black hole 
solutions have great achievements and provide a simplified model to find some conceptual issues in 
the context of black hole physics, quantum gravity, string and gauge theory, and also in the context of 
the AdS/CFT correspondence [50[feTl[521[551[Mll3?lf55] . In addition, critical behavior of the BTZ black 
hole with torsion has been considered in Ref. m- Also, Gregory-Laflamme instability of BTZ black 
hole in the context of massive gravity has been reported by Moon and Myung |58| . Thermodynamic 
equilibrium states of a typical static thin shell in a (2 + 1)—dimensional spacetime have been analyzed 
in [59] and thermodynamical properties of a BTZ black hole solution with an Horndeski source have 
been investigated in Ref. [60]. For more subjects related to BTZ solutions, we refer the reader to Refs. 

and references therein. 

Recently, the effective potential structure of the neutral BTZ black hole in gravity’s rainbow was 
studied in Ref. m- In this paper, we develop BTZ solutions to the case of existence of electrodynamics. 
We take into account the linear Maxwell field as well as some of the known nonlinear electrodynamics 
models in gravity’s rainbow. We also investigate thermodynamical behavior of the solutions and discuss 
thermal stability of them. 


2 BTZ and asymptotically BTZ Gravity’s Rainbow 


The Lagrangian of Einstein gravity with cosmological constant coupled to an electromagnetic field may 
be written as 

L tot =R-2A-C(F), ( 8 ) 

where R and A are, respectively, the Ricci scalar and the cosmological constant, and electromagnetic 
Lagrangian, C(F), is a function of Maxwell invariant T = F a bF ab , where the Faraday tensor is F a b = 
d[ a Aty and Ab is the gauge potential. Taking into account the gauge-gravity Lagrangian ([8]) and applying 
the variational method, we obtain 

G a b + Ag a b = T a b, (9) 

V a {C^F ab ) = 0, (10) 

where G a b is the Einstein tensor, Cjr = dC d ^ and 

T ab = \g^L{R) - 2 F^FfL?, (11) 


In this paper, in addition to the linear Maxwell electrodynamics (MED), we consider four classes of 
NED fields, namely BI nonlinear electrodynamics (BINED), Exponential form of nonlinear electrody¬ 
namics (ENED), Logarithmic form of nonlinear electrodynamics (LNED) and Power-Maxwell invariant 
form of electrodynamics (PNED) in which their Lagrangians are [711l75H76li77H751f791[5ni[5Tll821l831[5^1 
[85] 

-T, MED 


L{F) = 


4/3 2 (l - ^ 1 + , 

/3 2 (exp(-^r) - l) , 
-8/3 2 In (l + ^ 2 ) , 


(-JT 


BINED 

ENED 

LNED 

PNED 


( 12 ) 
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In these equations, /3 and a are called the nonlinearity parameters. It is notable that for /3 — > oo and 
a — > 1, all the mentioned nonlinear models reduce to Maxwell Lagrangian. 

Before we proceed, we would like to provide some reasonable motivations for considering these 
forms of NED. 

Although Maxwell theory is capable for describing different phenomena in classical electrodynamics, 
it is not flawless and has some fundamental problems [M115T1 . In order to solve these problems, we 
can generalize Maxwell theory to the case of NED [TSIl 77II7 8I[7^1[501I5I1[5^1I551IM1[551[5S1[?71I55] . Another 
strong motivation of considering NED, comes from developments in string/M-theory. It has been 
shown that the Born-Infeld [74] (Bl)-type theories are specific in the context of NED models, which 
are naturally acquired in the low-energy limit of heterotic string theory [59MMTH92] . In addition, 
it is easy to show that the weak field limit of BI type models is the same as calculation of the one- 
loop approximation of QED, which leads to quadratic Maxwell invariant in addition to the Maxwell 
Lagrangian [93] . 

Recently two BI -types of NED have been introduced, in which we can remove the divergency of 
the electric field of point-like charge near the origin. It was shown that Bi-types of NED enjoy special 
and important properties. To name a few, one can indicate the absence of shock waves, birefringence 
mm and an electric—magnetic duality 196] . It was shown that in a uniform electromagnetic field 
background for electrons, calculation of exact one-loop corrections leads to presence of logarithmic 
form of electromagnetic field strength Ezj • This correction term comes from the vacuum polarization 
effects. On the other hand, Bi-type theories have been employed for explaining the equation of state 
of radiation for inflation [881 . In addition, in the context of the holographic superconductors, it was 
shown that these NED models make significant effects on condensation, critical temperature of the 
superconductor and its energy gap [see [981199111001110111102] for more details]. It was shown that the 
exponential form of NED has stronger effect on the condensation formation and conductivity comparing 
to logarithmic form of NED [101111021 . 

Now let us present some motivations of considering PMI theory. The PMI theory is significantly 
richer than the Maxwell field, and in the special case (s = 1) it reduces to linear electrodynamics [521 
mmm- One of the most important properties of the PMI model in d-dimensions occurs for s = d /4 
where the PMI theory becomes conformally invariant and so the energy-momentum tensor will be 
traceless (the same as Maxwell theory in four-dimensions 179118011811 ]). In this case, one can obtain an 
inverse square law for the electric field of point-like charge in arbitrary dimensions. Furthermore, it 
was shown that there is an interesting relation between the black hole solutions of a class of pure F(R) 
gravity and those of conformally invariant Maxwell source in Einstein gravity 1103] . In the context 
of AdS/CFT correspondence, the effects of PMI source on strongly coupled dual gauge theory have 
been investigated [1041 (for more motivations of PMI theory, we refer the reader to [ 105] ]). Motivated 
by the recent results mentioned above, we consider the linear and nonlinear electrodynamics in three 
dimensional gravity’s rainbow. 

It was shown that adding the curvature to doubly special relativity (DSR) [T8] leads to an extension, 
called doubly general relativity (DGR) [35]. In DGR the spacetime geometry depends on the energy 
of the particle (E). Thus, one may regard a spacetime which is parameterized by the energy ratio e to 
form a rainbow of metrics. Three dimensional static metric in the framework of gravity’s rainbow is 
given by [35] 


dr 2 = — ds 2 


*( r ) ]f 2 _ 

P(e) g*(e) Uto 



(13) 


In this paper, we are looking for the black hole solutions with a radial electric field. Thus the 
nonzero components of the gauge potential is temporal component. Regarding the mentioned models 
of Eq. (fT91) and using the gauge potential A M = h(r)6®, we find that the electromagnetic field equation 
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m leads to the following differential equations 

rh" + h' = 0 , 


rh" + ti 


1 - 




1 + 


^ 2 /(e)g(e)ft' y 


MED 

= 0, BIN ED 

h" + h' = 0 , EWE 1 !? 


(14) 


rh' 1 


1 - 


1 + 

^ f(e)g(e)h' y 


+ 

b! = 0 , 


LNED 


r (2a - 1) h" + h' = 0, (a ^ 1) P1VEE 

where prime and double primes denote first and second derivative with respect to r, respectively. The 
solutions of the mentioned differential equations are 


h(r ) = < 


^ ln (f) > 

q\n[%(l + r)], 

'p(-¥) + -M J 


^ 2 r 2 (r-i) „ , 

qTWWF 2+« ln 


r(i+r) 

21 


( 2s ~ 1 )g 

2(1 — a) ) I 5 P D 

2(s-l)r 


MED 

BINED 

ENED 

LNED 

PNED 


(15) 


/ f2 (V)q2 (f)q2 

where E = w 1 + J //A — , 9 is an integration constant which is related to the electric charge of the 
black hole. In addition, L\y = LambertW ) which satisfies LambertW(x ) exp [Eam&erfIT(x)] 

o° 

x, 7 = 7 ( 0 ) ~ 0.57722 and the special function Ei (l,x) = f dz (for more details, see [1061 1. 

1 

According to the gauge potential ansatz with Eq. (1151) . one finds that the non-vanishing components 
of the electromagnetic field tensor can be written as 


Ftr — Ert. — < 


£ 
r 5 

MED 

Q 

BINED 

r r ’ 


fexp(-^), 

ENED 

«/W(s) ( r 1 ^’ 

LNED 

Q 

_JL_ j 

PNED 


(16) 


In order to have a sensible asymptotic structure for the electromagnetic field, F tr should vanish for 
large values of r. This condition is satisfied for MED, BINED, ENED and LNED branches. For PNED 
case, vanishing F tr for r — > 00 leads to a > 1/2. Thus in this paper, we restrict ourselves to a > 1/2. 
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Regarding the gravitational field equation for the mentioned models of NED, we find that the following 
metric function satisfies all components of the field equation ([5]), simultaneously 


where 


0(r) = < 


Ar 2 

^( r ) = ~~g(EY ~ m + 


-2g 2 / 2 ( e )ln(i), 

V/ 2 (e)ln ( r ^TT- i ) , 
q 2 P(e)E l ( l,Y) + |^y- 


g 2 f 2 (e)( 7-3+ln2) 
a 2 (A 


qPrf(c) 

g{e)VLW 


2g 2 / 2 (e)ln(^M£)£) 
4 q 2 f 2 {e) exp 


' —Lv 


f 2 {e)q 2 {- In 

6/3 2 r 2 r + 2} 


2“[/ 2 (e)9 2 (e)9 2 l“(2a-l) 2 
2(ct—l) 9 2 (e)r 2a_1 


MED 


BINED 


ENED 


LNED 


PNED 


and to is an integration constant which is related to the total mass (M) 

TO 


M = 


8/(e)' 


(17) 


(18) 


(19) 


It is notable that the ArnowittDeserMisner (ADM) mass of black hole can be obtained by using 
the behavior of the metric at large r. Straightforward calculations show that ADM mass of the black 
hole solutions is the same as Eq. 09- 

Using series expansion of metric functions for large values of r, we find that A-term will be domi¬ 
nated and therefore obtained solutions are asymptotically AdS with an effective cosmological constant, 
A e ff = g (E \2 ■ In addition, calculations show that the Ricci and the Kretschmann scalars for the metric 
(Thill are 


R = -g 2 (e)*" - 


2 g 2 (e)r 


R^pvR^ = g\e)E" 2 + 2 


g 2 (eW 


( 20 ) 

( 21 ) 


Inserting Eq. 03 into the Eqs. m and 03, we find that there is an essential singularity located 
at r = 0. Depending on the values of free parameters, one confirms that the mentioned singularity 
can be covered with an event horizon Hi- Applying the surface gravity interpretation on the event 
horizon, one finds the Hawking temperature as 


T = - 


Ato 


2 Trf{£)g(e) 


T. 


+ J 


( 22 ) 
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where 


-/(e)9(e)9 2 

27rr_|_ ’ 

MED 

r/3 2 (l -r+) 

’'7(e)ff(e) ’ 

BINED 

9/5(1 -L w+ ) 

P 2 r+ 

ENED 

2-KyJ Lw + 

4tt/( e)g(e) ’ 

2/3 2 r + [l+ln(^T±)-r + ] 

LNED 

*7(e)g(e) ’ 

-[2/ 2 ( E ) g 2 ( £ ) ? 2 ]“(2a-l) 

1 5 

PNED 


4?r/(e)a(e)r + “ 


(23) 


and r + = y 1 + ^^ fp{ e)q2 and L w+ = LambertW ( 4/ ^ r £ e)q )■ 

It is worthwhile to mention that in the absence of rainbow’s effects (/(e) = g(e) = 1), there is 
a critical value for the horizon radius (r+ c ), in which the temperature is positive for r+ > r+ c . The 
value of r +c depends on the choices of other parameters. In order to have a consistent limitation, we 
restrict /(e) and g(e) to smooth positive functions. It is evident that due to the presence of rainbow’s 
functions, the behavior of temperature and the values of r +c change, drastically. 


3 First law of thermodynamics and thermal stability 

In order to check the First law of thermodynamics, we should calculate entropy, electric potential and 
charge. Since we are working in Einstein gravity, it is allowed to use the area law to obtain the entropy 
[107U1081I1O9] . The entropy of black holes is equal to one-quarter of the horizon area 


7 rr + 

Me)' 


(24) 


Besides, we can obtain the electric charge by using the Gauss law. By calculation of the flux of the 
electric field at infinity, one can find the electric charge as 


Q = < 


g/(e) 

2 > 


MED 

BINED 

ENED 

LNED 


(25) 


2°/q 
. 9(e) 


[qf{e)g{e)] 2a -\ 


PNED 


where for a = 1, the electric charge of PNED branch is compatible with other branches. The electric 
potential measured at the reference with respect to the horizon, is defined by 


^ |reference r+ ) 


( 26 ) 























where L> vanishes at the reference and x^ is the null generator of the horizon, yielding 


L> = 


-9ln[^(l + r + )], 
l-exp(— 


In ( ^ /( g g(e)g ) - 


& 


lP{s)g' 2 {e) 


— gin 


(i+r + )r + 

21 


" (2 “ 2 I 1 I) , (a / 1) 


MED 

BINED 

ENED 

LNED 

PNED 


(27) 


2(a-l)r + 2 ' 


After calculating all of the conserved and thermodynamic quantities of the solutions, we are in a 
position to examine the first law of thermodynamics. To do this, we obtain the total mass M as a 
function of the extensive quantities Q and S. Using the expression for the entropy, the electric charge 
and the mass given in Eqs. m , dm (Ei and El. and the fact that f(r = r+) = 0, one can obtain 
a Smarr-type formula as 

AQ2 y 

= (28) 


2 7r 2 /(e) /(A)’ 


where 


E = 


Q 2 to (*§£), 




Ql 

2 

Ql 

2 


Ei 

In 


i ( 1 >^)+'r-3 

(8 QV(£)\ 

[ P 2 l 2 ) " r 4 tt 2 

/3QS(l-2L' w ) 


MED 


BINED 


ENED 


(29) 


27Ta/L 


'w 


21 


(3—r' 2 )in(a±£l) (2 _ r>) 

( 1 -r' 2 ) + (l+r') -1- 


i„(M4) 


LNED 


PNED 


(2ct-l) 2 2 “ _1 

2 —2a 5 

4(a-l)(f)5^T 

and r' = 1 + and = LambertW { ^^2 32 Taking into account Eq. 

the temperature and electric potential as intensive parameters 


\dSj Q ’ VWs 

which are the same as those calculated in Eqs. (HE1) and ESI) (fo: 
that obtained conserved and thermodynamic quantities satisfy 1 


we can calculate 

(30) 


>r all branches). Therefore 




dM = TdS + <L>dQ. 


(31) 
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In what follows, we investigate thermal stability. In principle, thermal stability of a system with 
respect to small variations of the thermodynamic coordinates is usually performed by investigating the 
behavior of the entropy S(M, Q) around the equilibrium. 

The local stability criteria in any ensemble requires that S(M,Q) be a smooth convex function 
of the extensive variables or its Legendre transformation should be a smooth concave function of the 
intensive thermodynamic coordinates. Thermal stability can also be studied by the behavior of the total 
mass M(S, Q) which should be a smooth convex function of its extensive variables. Therefore, the local 
stability can in principle be carried out by obtaining the determinant of the Hessian matrix of M(S , Q ) 
versus its extensive coordinates X,, H ^. x = [d 2 M/dXidXj] (see Refs. [1101111111112] for more details). 
In our case the mass M is a function of the entropy S and the charge Q. The number of thermodynamic 
variables depends on the ensemble that is used. In what follows, we calculate the heat capacity of 
(nonlinearly) charged black hole solutions to investigate thermal stability through canonical ensemble. 
We should regard the electric charge as a fixed parameter and calculate Cq = T (§§)q- The positivity 


of the heat capacity guarantees the local thermal stability. Since we investigate the physical black hole 

-l 


solutions (positive temperature), it is sufficient to examine the positivity of 


<dS\ 
\dT ) Q 


( d 2 M \ 


for 


thermal stability. Calculations show that 


d 2 M\ —A 


+ 


1 /(e)g (e) 

7 r 2 r'^_ ' 


2/3 2 (r + -i) 

27T 2 /(e:) 

4/3 2 ln(^) 
ir 2 f(s) > 


2°l9/(£)g(e)l 2 ° 

2 TrS/OOrf^ 


MED 

BINED 

ENED , 

LNED 

PNED 


(32) 


L W+ { Ft 2 m\ 

It is clear to find that T + and e~~ are greater than one and therefore I 1 is positive definite 

for obtained AdS solutions. Regarding positive values for f(e ) and g{e) and taking into account a 
positive temperature for the physical black hole solutions ( for r+ > r+ c ), we can conclude that the 
mentioned black holes enjoy thermal stability. In comparison with the asymptotically AdS black holes 
of Einstein gravity, which have a small unstable phase, the stability phase structure of the Einstein- 
nonlinear) Maxwell black holes in gravity’s rainbow shows that the electric charge help us to obtain 
stable solutions and rainbow functions do not make anomaly. 


4 Note on the integration constants 

Regarding positive real rainbow functions with the functional forms of r and L\y (after Eq. (1151) 1. one 
finds that there is no restriction on q and m. Besides, one can follow the method of Refs. mm 
HT51 with a specific model of rainbow function to obtain possible mass/charge limitations. Since in this 
paper we have discussed general form of rainbow functions, we release the mentioned method. 

Taking into account Eqs. m and m, one can find that q and m are integration constant. Since 
the variable of integrations is r coordinate, we can replace q and m with modified constants such 
as qf a (e)g b (e ) and mf c {e)g d {e ), respectively (a, b, c and d are arbitrary real valued numbers). It is 
worthwhile to mention that although in new situation, all conserved and thermodynamic quantities will 
be modified, the results for e —> 0 (/(e) = g(e) = 1) do not change. Another comment which should be 
discussed is the first law and thermal stability. Regarding new situation (applying q —> qf a {e)g b {e) and 
to —> mf c {e)g d {£) in previous sections), it is easy to show that the first law is satisfied. In addition, the 
criteria of thermal stability do not change, globally. These behaviors are expected, since we applied a 
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global transformation on the constants. Considering new situation, one can obtain four free parameters 
a, b, c and d, and their manipulations with special choices may lead to interesting results. 


5 Closing Remarks 

In this paper, we took into account the gravity’s rainbow in three dimensions in the presence of linear 
and nonlinear models of electrodynamics. The motivation of considering gravity’s rainbow comes from 
the quantum gravity view point, while the reason of considering nonlinear electrodynamics is arisen 
from LQG as well as low effective limit of heterotic sting theory. We found that although the rainbow’s 
functions affected the electromagnetic field behavior of nonlinear Bi-type models, Maxwell and PNED 
models of electromagnetic field did not depend on the rainbow’s functions. We obtained metric functions 
of each models and found that these solutions can be interpreted as black holes. 

Then, we calculated conserved and thermodynamic quantities and found that these solutions may 
depend on the choice of rainbow’s functions. We also considered finite mass as a function of extensive 
parameters and checked the first law of thermodynamics. In addition, we analyzed thermal stability 
by using of the heat capacity. We found that obtained physical solutions (T > 0) enjoyed thermal 
stability in the canonical ensemble. 

It is worthwhile to generalize these three dimensional solutions to four and higher dimensional cases 
with various horizon topologies. It is also interesting to investigate the effects of rainbow’s functions 
on dilaton gravity as well as higher curvature theory. These cases will be addressed elsewhere. 
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